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We investigate the electromagnetic radiation released during the high energy collision of a charged 
point particle with a four-dimensional Schwarzschild black hole. We show that the spectra is flat, and 
well described by a classical calculation. We also compare the total electromagnetic and gravitational 
energies emitted, and find that the former is supressed in relation to the latter for very high energies. 
These results could apply to the astrophysical world in the case charged stars and small charged black 
holes are out there colliding into large black holes, and to a very high energy collision experiment in 
a four- dimensional world. In this latter scenario the calculation is to be used for the moments just 
after the black hole formation, when the collision of charged debris with the newly formed black 
hole is certainly expected. Since the calculation is four-dimensional, it does not directly apply to 
Tev-scale gravity black holes, as these inhabit a world of six to eleven dimensions, although our 
results should qualitatively hold when extrapolated with some care to higher dimensions. 
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I. INTRODUCTION 

In a previous paper we have investigated the gravi- 
tational radiation emitted by a relativistic particle in- 
falling into a four-dimensional Schwarzschild black hole 
0, and afterwards extended it to the infall in a Kerr 
metric spacetime Q- The interest of these studies lies 
in the fact they can describe the gravitational radiation 
emitted in several collision phenomena, such as the col- 
lision between small and massive black holes, between 
stars and massive black holes, or between cosmic rays 
and small black holes, to name a few. The possibility of 
detecting gravitational radiation by the several operating 
antennae is now real, and a closer understanding of these 
collisional processes is accordingly important. Of course, 
to fully understand collisional processes, one has in prin- 
ciple to go beyond the infall of a particle and consider 
the collision between two black holes with comparable 
masses, in which case one has to take into account the full 
non-linearity and strong field regime of Einstein's equa- 
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tions. Without resorting to numerical computations, but 
instead using topological arguments, Hawking was able 
to put upper limits on the gravitational radiation emitted 
in a black hole - black hole collsion in four dimensions. 
This calculation was then refined by D'Eath and Payne 
3]. Their calculation made use of the fact that, if one 
boosts the Schwarzschild metric to high velocities, then 
it approaches the Aichelburg-Sexl Q metric, which is a 
shock wave spacetime describing the gravitational field 
of a massless particle. It was found that the total effi- 
ciency of the shock would be ~ 16% (see Yoshino and 
Nambu 5] for the analysis in higher dimensions). In 
the inclusion of a second term in the news function brings 
about a decrease in the efficiency from 25% (using only 
the first term) to 16%, so one is entitled to ask what de- 
gree of confidence do these results and techniques offer. 
For instance, one might suspect the third term to lower 
even more the efficiency. One needs to have other means 
of computing this process, since D'Eath and Payne's for- 
malism is hard to pursue, due to its complexity. We have 
recently 0, Q proposed the "point particle paradigm" , 
which may well be a good candidate. Handling high en- 
ergy collisions of point particles with black holes is rather 
simple, since the point particle may be looked at as a 
perturbation in the black hole spacetime, and the full 
machinery to handle black hole perturbations is well de- 
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veloped 6]. The interesting thing about the high energy 
collision of point particles with black holes is that in the 
limit that the mass /z of the particle goes to the mass M 
of the black hole, the result for the total radiated grav- 
itational energy agrees well with other predictions and 
never violates the area theorem 0,0- Moreover, one can 
use this method to treat the collision of rotating holes, 
which does not seem feasible using previous techniques. 

Up to now, most studies have concentrated on the 
emission of gravitational radiation in the collision pro- 
cess, but if at least one of the particles is charged, elec- 
tromagnetic radiation is also expected. Astrophysically 
this is perhaps less relevant for the time being, since the 
particles one can consider, such as stars or small black 
holes, are in principle discharged, although it is possi- 
ble the latter ones can carry some charge (see, e.g., Q). 
On the other hand, from the point of view of elementary 
black holes such a study is certainly of interest: elemen- 
tary black holes and particles can easily carry electromag- 
netic charge, and the collision between such objects can 
be put to operate, either by a machine in some far future 
that can collide particles at center of mass energies of 
the order of the four-dimensional Planck mass, 10 16 TeV, 
or if there is any fundament in the TeV-scale gravity, by 
the near future generation accelerators such as the LHC. 
The latter case has been given a lot of consideration. The 
TeV-scale gravity requires extra large dimensions (of 
the order of sub-millimeters or smaller) in order to lower 
the higher-dimensional Planck mass to energies of the 
order of TeVs. In these TeV-scale gravity models, the 
gravitons are free to propagate in the higher dimensional 
spacetime, whereas, due to experimental constraints, the 
standard model fields live on a 3-brane, our universe. In 
this scenario, a particle collider with a center of mass en- 
ergy of the order of TeVs can copiously produce higher 
dimensional small spherical type black holes, with radius 
of the order of fermis or less, in a space with large extra 
dimensions Q, or even black branes if some of the ex- 
tra dimensions are large and other small (see [T5 ] for a 
review). Moreover, in the collision process the gravitons 
escape to the extra dimensions and are therefore much 
harder, harder than usual, to detect. On the other hand, 
all the electromagnetic radiation emitted in such a colli- 
sion can be detected. 

Thus it is important to know the electromagnetic spec- 
trum and the quantity of electromagnetic energy radi- 
ated in such an encounter. In this paper we extend 
the previous calculations into to the electromagnetic 
window and find the electromagnetic radiation emit- 
ted by a highly relativistic electrically charged parti- 
cle infalling into a four-dimensional Schwarzschild black 
hole. This is a perturbation calculation: we suppose a 
small charged particle falling into a Schwarzschild black 
hole. Some former works that have dealt with the phe- 
nomenon of electromagnetic radiation from a charged 
particle falling from infinity into a Schwarzschild black 
hole are 

[H Ellll For example Ruffini [H first cal- 
culated the electromagnetic energy spectra and total ra- 



diated energy |12(, for particles with low Lorentz factors, 
7 << 1. Here we go into high 7s. We use the point parti- 
cle approximation and consider a charged point particle 
colliding head-on at high velocity with a Schwarzschild 
black hole. We show that the spectra is flat, and we also 
compare the total electromagnetic and gravitational en- 
ergies emitted, and find that the former is supressed in 
relation to the latter for very high energies. The numer- 
ical results extracted by us are in very good agreement 
with Ruffini's results [IJ. We shall also see that there 
is a classical calculation for this process that agrees ex- 
tremely well with our numerical results. 

Two comments are in order: (i) Our calculation is for 
the collision of a small particle with a black hole in a 
four-dimensional world. Therefore, in principle, it could 
apply to the astrophysical world in the case charged stars 
and charged black holes are out there, and to a very 
high energy collision experiment in a four-dimensional 
world (without extra dimensions) and with a fundamen- 
tal Planck mass of 10 16 TeV. In this latter scenario the 
calculation is to be used for the moments just after the 
black hole formation, when the collision of charged debris 
with the newly formed black hole is certainly expected. 
Furthermore, the calculation can be extended to a black 
hole - black hole collision as we have already argued 0,0, 
and, in addition, it can give clues, although it does not 
aplly directly to the usual collision process in a collider, 
i.e., the collision between one non-charged particle with 
strong gravitational field (not necessarily a black hole) 
and one charged particle. These results may, however, 
serve as a model to the electromagnetic radiation emit- 
ted in the initial phase of a newly formed black hole, the 
stage in which the black hole sheds its hair by emitting 
gauge radiation, such as electromagnetic radiation, (ii) 
Our calculation is four-dimensional. It does not apply to 
Tev-scale gravity black holes, since these inhabit a world 
of six to eleven dimensions, with some of the dimensions 
being large, others perhaps being small. However, al- 
though there are certainly some differences, qualitatively 
our results should hold when extrapolated with some care 
to higher dimensions. 

II. BASIC FORMALISM 

In this study, we assume the charged particle and 
the emitted radiation to be a small perturbation on the 
Schwarzschild spacetime, whose line element is given by 

ds 2 = -e v dt 2 + e x dr 2 + r 2 (d9 2 + sin 2 6 dip 2 ) , (1) 

where e v — e _A = 1 — 2M/r, and M stands for the mass 
of the background spacetime. Because a complete pre- 
scription to deal with this problem was given by Zcrilli 
[Tl| , here we briefly show the essential parts of the for- 
malism. In virtue of the static property and spherical 
symmetry of the background spacetime, perturbations 
can be decomposed by using vector harmonics for the 
angular variables 9 and tp, and the Fourier components 
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for the time variable t. After accomplishing a separation 
of variables, perturbations are therefore characterized by 
the harmonic indices I and m, parity, and frequency u>. 
For the present case, where a particle falls straight into 
a black hole along the z-axis or the 9 — line, no axial 
parity perturbations are excited due to the symmetry of 
the motion. Thus, only the polar parity perturbations 
are considered in this study. According to Zerilli [Tl| . 
our master equation for determining the electromagnetic 
radiation is given by a single wave equation, 



d 2 f lm (u),r) 
dri 



1(1 + 1) 



fl m (cj,r) = e"Si m , (2) 



where r* is the so-called tortoise coordinate, defined by 
dr/dr* = e^-^/ 2 , or r* = r + 2Mlog(r/2M - 1) in 
our case, and Si m is the source term determined by the 
charge and motion of the particle. In the case of the 
radially falling particle, Si m is generally given by 



1 e iwT(r) 
~ 2( l\l 1 + « 9 S m0 ■ 



(3) 



where <5y means the Kronecker delta, and q denotes the 
charge of the particle. Here, the function T(r) is the co- 
ordinate time of the particle parameterized by the radial 
position r of the particle. Since the charged particle is 
treated as a perturbation of the background spacetime, 
the particle traces a radial geodesic of the Schwarzschild 
geometry. Thus, the function T(r) is determined by (see 
e.g., Chandrasekhar [T^ l 



dTjr) 
dr 



(4) 



where 7 is an integral of motion, given by 7 = (1 — 
v 2 ^) -1 / 2 , where Woo is the radial velocity of the particle 
at spatial infinity. 

Once solutions fi m (u,r) of equation are calcu- 
lated, time dependent functions fi m (t, r) can be obtained 
through the inverse Fourier transformation, 



flm{t,r) = 



t fi m (u,r)duj . 



(5) 



From the functions fi m (t,r), all the components of the 
field strength F a p can be derived except at the position 
of the charged particle. For example, the 1 — 6 component 
F r g is given by 



F re (t, r, 6,<p) = e- v Y, r )-£r( 9 > <P) • 



I . m 



(6) 



In order to obtain a unique solution of equation J2J), 
two boundary conditions must be specified, and physi- 
cally acceptable conditions are purely outgoing wave at 
spatial infinity and purely incoming wave at the horizon, 
which are due to equation JSJ) given by 



fim(r) 



Irn 

A^(co)e^ 



cLS T* 



-00 , 



00 . 



(7) 



where A°^(uj) and are not dependent on r. Since 

we are interested only in solutions at spatial infinity, what 
we have to do is to obtain Af^(aj)'s as functions of u for 
our purpose of this study. Adapting a standard Green's 
function technique, we can write as the integral, 

given by 



AZ'H = 



1 



2iujCi(uj) J 2M 



fh,iSi m dr , 



(8) 



where /l,z is a homogeneous solution of equation J5J) sat- 
isfying a boundary condition given by 



/l.i 



e * as r* — > — 00 , 

Bi(uj) e lujr - + Ci{uj) e~ wr * as r* -> 00 . 



The energy spectrum at spatial infinity is given by 



(9) 



dE v ^ dEi x - 



l{l + l) 
2tt 



|A° ut H| 2 for c>0, 
( 10 ) 

In order to obtain numerical values of A°^(uj), we start 
the integration of /lj and J^ M fh,iSi m dr' at r = 2M(1 + 
10 -6 ) by using a Runge-Kutta method. Those functions 
are then integrated out to large values of r. The inte- 
gration is stopped if the absolute value of J 2M fajSimdr 1 
converges within the required numerical accuracy, and 
we simultaneously match faj with an asymptotic solu- 
tion satisfying condition 10 at spatial infinity, which is 
given in [lj|, to obtain a value of Ci(uj). 



III. NUMERICAL RESULTS 

Following the numerical procedure just outlined, we 
have computed the energy spectrum for several values of 
7. In Figure 1 we show a typical result of the energy 
spectrum, here for 7 = 2 and for the four lowest radiat- 
ablc multipolcs. In Fig. 2 we show similar results but 
for a particle with 7 ^ 10. A general feature of the 
energy spectrum for high energy collisions is that it is 
flat up to some critical frequency, after which it rapidly 
(exponentially) decreases to zero. This was also verified 
for the gravitational energy spectrum resulting from the 
high-energy collision of a black hole with a point par- 
ticle 0, y|- This critical frequency is given, in a very 
good approximation, by the frequency of the fundamen- 
tal quasinormal mode associated with the corresponding 
harmonic index I. The corresponding quasinormal fre- 
quencies are also indicated by a filled circle in Figures 1 
and 2. The total energy AE radiated away is given by 



AE 



E 

1=1 



AE, 



(11) 



where AEi is the total radiated energy associated with 
Z-pole electromagnetic radiation, which is given by the 
area under the curve of the energy spectrum dEi / du (see 
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FIG. 1: The electromagnetic energy spectra for the four low- 
est radiatable multipoles, for a 7 = 2 particle falling from 
infinity into a Schwarzschild black hole. The filled circles on 
the curves of the energy spectra indicate the frequency of 
the fundamental quasinormal mode associated with the cor- 
responding harmonic index I. 
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FIG. 2: The electromagnetic energy spectra for the four low- 
est radiatable multipoles, for a 7 ^> 10 particle falling from 
infinity into a Schwarzschild black hole. The filled circles on 
the curves of the energy spectra indicate the frequency of 
the fundamental quasinormal mode associated with the cor- 
responding harmonic index I. 



TABLE I: The total radiated electromagnetic energy as func- 
tion of 7 . 



7 



AEMq~ 2 
2.14 x 10" 2 
1.14 x 10" 1 
2.23 x KT 1 
3.34 x KT 1 
4.46 x KT 1 



Figure 1). For the low 7 cases, the convergence of the 
infinite series (|11J) is relatively good in the sense that a 
finite summation of several AE^s is sufficient to obtain 
a value of AE with sufficient accuracy. However, the 
convergence of the total radiated energy gets worse with 
increasing 7. This contrasts with the analogous calcula- 
tion of the gravitational energy radiated in high energy 
collisions Q where a rescaling was possible to show that 
AE cx 7 2 . Here we shall see in the next section that the 
difficulty for the electromagnetic case lies in the fact that 
AE oc log 7, for high 7's. For a very fast particle, 7 ~ 00, 
it is quite difficult to numerically obtain the total radi- 
ated energy, since all the Z-pole energy spectra dEi/dw 
make a substantial contribution to the total radiated en- 
ergy. In this study, we calculated the total energy for 
several values of 7 with 7 < 10. Some of these values are 
shown in Table 1. 

Previous works concerned with the electromagnetic ra- 
diation from a charged particle falling from infinity into 
a Schwarzschild black hole include for example Ruffini 
p^ |. Ruffini first calculated the electromagnetic energy 
spectra and total radiated energy f° r l° w 7 particles. 
The numerical results extracted bv_ us are in very good 
agreement with Ruffini's results 



IV. THE CLASSICAL CALCULATION 

There is a classical calculation in electromagnetism 
that suggests itself as a model to use in the high energy 
collision of charged particles: the radiation emitted when 
a charge, with constant velocity Voq is suddenly deceler- 
ated to v = 0. The deceleration is idealized as taking zero 
seconds. This model calculation has been applied with 
great success to, for example, beta decay [Tj, Ha. The 

is~|l!G3 



result for the energy spectrum per solid angle is 



1 1 



sin 



4tt 2 (1 



\dujdQ J 
or, integrating over solid angle 



dE 
du) 



2 r 



flog f 

Woo V 1 



- 2 



(12) 



(13) 



This formula has been tested with great accuracy exper- 
imentally. To get the total energy, one has to integrate 
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(|13fl over frequencies, and a naive procedure would lead 
to infinities in the total energy. To obtain reasonable re- 
sults, one has to impose a cutoff uj c in the frequency de- 
pending on the particular problem under consideration, 
and one obtains 



1 , fl + V 
l0g 1 



- 2 



(14) 



How well do these classical formulas fit into our numer- 
ical results? Very well indeed. To allow a more direct 
comparison, we shall first decompose the energy spec- 
trum Ijl3|) into spherical harmonics, i.e., into the energy 
spectrum associated with each Z-pole radiation as follows: 



dE 
du) 



:'la.ss 
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(15) 
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When we consider the limit of Voo — > 1 , as shown in equa- 
tion (|13|) . the total energy spectrum (^fj) class diverges. 
However, the energy spectrum due to Z-pole radiation 
(4^") class converges to a finite value, given by 



FIG. 3: The cutoff frequency, given as a function of 7. The 
filled circles on the curve of the cutoff frequency indicate the 
cutoff frequency calculated in this study. The linear function 
uj c = 0.224 + 0.0598 7 is also shown. 



dEA 
duj ) 



2 2/ + 1 

dE 



(16) 



This yields for example , = 0.477g : 
0.265g 2 as 7 — > 00. This may be compared with the nu- 
merical values in Fig. 2. The numerical agreement is 
excellent not only for these lower multipoles, but also for 
higher ones, and one has therefore established numeri- 
cally that for high energy collisions the classical result 
(|13fl for the energy spectrum is valid. Probably more in- 
teresting is the total energy radiated. Since one already 
knows that the classical result for 4^ is correct, let us 
now try to predict the cutoff frequency by identifying 
AE in equation i|TT|) with ljT4"|> . Making use of this defi- 
nition of lo c and numerical values of AE, we evaluate the 
cutoff frequencies cu c in (|14|l for several values of 7, and 
list the values of oj c in Table II. In Figure 3, we also show 
the cutoff frequencies as a function of 7. It is found in 
this figure that for relatively high vales of 7, u> c increase 
almost linearly with 7, and a good fit to our numerical 
values is given by 



TABLE II: The cutoff frequency ui c as a function of 7 



LU C M 



OO 

3.44 x 10" 1 
4.03 x 10 _1 
4.64 x 10 _1 
5.23 x 10 _1 



or considering that Voo ~ 1 we have also 
2o 2 

AE=-j- [log 2 7 - 1] (0.224 + 0.05987) , (19) 



V. COMPARING THE AMOUNT OF 
GRAVITATIONAL AND ELECTROMAGNETIC 
ENERGY RADIATED 



Moj c = 0.224 + 0.0598 7. 



(17) 



In Figure 3, we have shown this linear function, and we 
can confirm that this linear function is in good agreement 
with our numerical values for high 7. This means that a 
good approximation to the total electromagnetic energy 
radiated is 



AE 



g 

Mtt 



— log 

Woo 



1 



(0.224 + 0.05987), 

(18) 



Previous studies on the high energy collision of point 
particles with black holes led Q, Q to the following con- 
clusion: if a high energy point particle of mass /i collides 
head on with a mass M non-rotating black hole, the total 
amount of gravitational energy radiated is 



AE g = 0.26^- 



(20) 



where 7 is again the Lorentz factor for the point particle. 
Furthermore, it was also found that the spectra is flat and 
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also well described by a classical calculation (here clas- 
sical means again working on a flat background). Now, 
for high 7's (JTHJ is 

^-7log2 7 , 7^00. (21) 

So, we get for the ratio electromagnetic energy over grav- 
itational energy, 



0.146(— ) , 7 — > 00 . 



M 7 



(22) 



This means that electromagnetic energy is supressed in 
relation to gravitational energy, for very high energies. 



VI. CONCLUSIONS 

We have computed the electromagnetic spectrum and 
total energy radiated during the high energy collision of 
a charged point particle with a Schwarzschild black hole. 
Our results show that the classical "instantaneous colli- 
sion" calculation gives very good results, in accordance 
with the full numerical ones. We have dealt only with 
zero impact (i.e., head-on) collisions, but these results 
are easily generalized to non head-on collisions. For ex- 
ample, one expects that the classical results [l5llT6| hold 
for such cases, and therefore that the total energy de- 
creases as one increases the impact parameter. We stress 
that the results presented here are valid for any particle 
charge q, as long as the total effective stress energy of 



the particle is small compared to the total energy con- 
tent of the black hole. Previous works Q, have shown 
that the high energy collision between two black holes of 
equal mass may be well studied through the collision of a 
point particle with a black hole, and then taking the limit 
of equal mass, although this is formally not allowed. It 
is our belief that the same may be done here. One needs 
however some other method to attack this problem, to 
confirm or disprove this claim. The investigation car- 
ried here can also be carried over to higher dimensions, a 
case which is of more direct interest to TeV-scale gravity 
scenarios. Since Standard Model fields inhabit a four- 
dimensional brane, what one would need to generalize 
this construction would be to modify the induced four- 
dimensional metric as was done for example in |17| . 
The generalization is straightforward. 
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